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Recently, th? 




scattering-marri formalism a^J^jfe assu^gp_tion ^quantug-cohergctt trajpjsj^rt. liigcpQty;-ast^'|^p this,^^ paper Q 
this result was obtainea using quasiclassical kinetic equation under the assumption of completely incoherent local 
scattering. Although all the impurity-scattering ey^^jWem^^nsidcTCd to be independent, at low frequencies the 
current-conservation law resulted in averaging oin^jnin^nicvffmrit^over the contact volume, thus establishing a 
sort of correlation in electron transport. In a sense, this appears to be equivalent to quantum-coherent scattering, 
in which case the whole impurity system is considered as a single scatterer. In both approaches, the low-frequency 
shot noise appears to be independent of contact geometry [Q, ||. However, this should not be the case for 
the finite-frequency noise because it should be affected by long-range Coulomb interactions. In particular, it should 
depend on the possibility for the charge to pile up in the contact, i.e., on its external capacity. The problem of 
frequency-dependent noise in diffusive mesoscopic contacts was addressed by Altshuler, Levitov, and Yakovets || 
for the case of quantum-coherent transport, where the frequency dependence of noise was determined by fermionic 
correlations, but no account of long-range Coulomb interactions was taken in this paper. 

In the present paper, we consider the effects of contact geometry on the frequency dependence of shot noise within 
the semiclassical incoherent-scattering approach. We consider the case where all its dimensions are much larger than 
the screening length Ao. The contact of length L is either a cylinder of circular section with a diameter 2ro or a 
plane-parallel layer of thickness oIq consisting of a metal with a high impurity content (see Fig. |l|). The electrodes 
are of the same section, yet the resistivity of their material is negligible. The contact is embedded in a perfectly 
conducting grounded medium, which is separated from its surface by a thin insulating film of thickness 5q and the 
dielectric constant £d- As will be shown below, this particular choise of contact geometry allows us to avoid solving 
Poisson equation in the surrounding medium and reduces the effects of environmental screening to frequency-dependent 
boundary conditions. 

The external circuit is assumed to have a large grounding capacity, which allows accumulation of the charge in it. 

Our consideration is based on the Boltzmann - Langevin approach first proposed in H . The nonuniform extraneous 
currents caused by randomness of electron-impurity scattering result in local charge-density fluctuations in the bulk 
of the contact. These fluctuations are effectively screened by the surface charge induced at the outer surface of the 
contact and in the surrounding medium (environmental screening) and at the contact - electrode interfaces (electrode 
screening). As will be shown below, the finite-frequency noise essentially depends on the dominating type of screening. 
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shot\noise in mesoscopic contacts became a subject of extensive study. One of the principal results 



was that in (j^tacts Pith's -strong elastic_scattering,_the low-frequency, shot- noise is-l/3-ofthc-fulL noise of classical 
Poisson process. This result was obtained almost simultaneously by different authors using different methods and, more 
importantly, difl erent physical^assumptions. Beenakker and Biittikcr Q obtained this result using the m ultichannel 



II. BASIC EQUATIONS 

In the Boltzmann-Langevin approach, the long-range Coulomb interaction is taken into account by fluctuations of 
charge density Sp and self-consistent fluctuations of electrical field <$E. Consider the case of strong and purely elastic 
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scattering. The Boltzmann-Langevin equation for fluctuations reads 



d_ v _9 gEv _9 
dt dr de 



5f + 51 = -eSEv^- + SJ ext , (1) 
oe 



where SJ ext is the random extraneous flux. The correlation function of these fluxes is given by the expression [Q 

(6J ext (p,r,t)6J ext (p',r',l/)) = S(r - r')S(t - t') 

(V E^(P'P + «i)/(p + <0[i - /(p)] + w(p + q. p)/(p)[i - /(p + q)]} 

-{W(p', P)/(P)[1 - f(p')] + W(p, p')/( - /(P)]}) , (2) 



where W^(p,p') is the probability of scattering from state p' to state p. The fluctuation of electrical field S~E in the 
right hand side of Eqn. (|l|) is determined from the Maxwell equation 

V SE = 4tt 8 p. (3) 

The fluctuations of the charge and current density are given by the expressions 

Sp(r,t) = eJ d 3 P Sf(p,r,t), <5j(r, t) = e J d 3 p v<5/(p, r, t). (4) 

Now we proceed to the hydrodynamic approach and obtain a closed set of equations for macroscopic quantities <5j and 
dp. As the impurity scattering is strong, we can split the fluctuation of distribution function into the parts symmetric 
and antisymmetric in the momentum space. Then we separate the antisymmetric part of Eqn. (Q) from its symmetric 
part. Note that the extraneous flux contains only the antisymmetric part because the electron- impurity scattering 
does not change the total number of electrons at a given point with a given energy. Multiply the antisymmetric part 
of (|l|) by ev, integrate it with respect to d 3 p, and then multiply both its parts by the elastic scattering time, r. If the 
characteristic times considered are much larger than r, one obtains 

Sj = -D-^-Sp + aSE + 5j ext , (5) 
Or 

where D = v 2 t/3 is the diffusion coefficient, a — e 2 NpD is the conductivity of metal, and <5j ext = er J d 3 pv SJ ext . 
Integrating the symmetric part of Eqn. ([!]) with respect to momentum, one obtains just the current-conservation law 

^Sp + WSj = 0. (6) 

Applying the operator V to both parts of Eqn. (g) and making use of Eqns. (g) and §), one obtains a closed 
equation for fluctuations of charge density in the form 

4- - DV 2 + 4iro) 5p = -V Sj ext . (7) 
at J 

In the left-hand side of this equation, the second term describes diffusion of electrons, and the third term describes 
the Coulomb screening of fluctuations. In principle, Eqn. (f7|) may be solved for each particular distribution of <5j ext , 
and then £E and Sj may be determined from Eqns. (Q) and (||), respectively. In the static limit, Eqn. (Q) describes 
the screening of an extraneous charge with the standard length Ao given by 

,_n 47TCT , n, T 

A 2 = — =4ne 2 N F . 

To complete the derivation, we must obtain the correlation function of extraneous currents <5j ext . As we are restricted 
to the case of strong impurity scattering, the distribution function may be considered as isotropic in the momentum 
space and dependent only on the coordinate r and energy e. Multiply Eqn. (||) by erv a and erv'p, where a and f3 
label vector components, and integrate it with respect to d 3 p and d 3 p' . As a result, one obtains the spectral density 
of extraneous currents in the form 

(6f a xt (r)Sj^(r')) u =4aS a0 S(r-r') J de /(e, r)[l - f(e, r)]. (8) 
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Because of smallness of Ao , the relationship between the extraneous currents and fluctuations of charge density in 
the bulk of the sample may be considered as local. Taking the Fourier transform of Eqn. ([?]) with respect to time, 
integrating it over the space, and making use of the Gauss theorem, one obtains: 



Sp = -(-iw + 4na) 1 V8'f 



(9) 



Note that the quasineutrality condition does not hold for fluctuations. Introduce the fluctuating potential 86 that 
satisfies the Poisson equation 



-AirSp. 



(10) 



Consider the boundary conditions for 86 at the outer insulated surface of the contact. The normal derivatives of < 
in the dielectric layer and inside the metal are related by the expression 



85(f) 
dn 



d5(f) 
dn 



= — 4ir5(Te 



(11) 



where 5a s is fluctuating surface charge density induced by the extraneous currents. On the other hand, this charge 
density satisfies the charge-balance equation 



j5<j x 



88(f) 
dn 



(12) 



As the thickness of dielectric layer is much smaller than the size of the contact, the electric field across it may be 
considered uniform so that 85(f>/dn\ d = —Sq 1 5(f>\ s . With this condition, Eqns. ( |TT1 ) and ( |T^ ) give the boundary 
condition for 56 in the form 



-iojEd5 1 86 + (—iuj + 47T(t) 



886 
dn 



(13) 



It is easily seen that at lo — 0, Eqn. (|i3| ) takes the form 856/8n\ s — 0, while at ui — > oo, it takes the form 86\ s = 0. 

As the voltage drop across the contact is held constant, fluctuations of potential are zero at the contact-electrode 
interfaces: 



0. 



(14) 



As the electrodes are perfect conductors, 86/ dn = inside them. Equation ( pT[ ) holds for contact-electrode interfaces, 
but the charge-balance equation takes now the form 



iuo8a s 



dn 



5jr 



(15) 



where 5j n is the fluctuation of current flowing into the electrodes from the contact. From Eqns. (|llj) and ([L5|), it 
follows that the density of outgoing current is given by 



Sjn = 



-in 



dn 



(16) 



From the standpoint of average current, the problem is purely one-dimensional, so the average distribution function 
f(e,x) obeys the one-dimensional diffusion equation, its boundary values being zero-temperature Fermi distribution 
functions shifted in energy by eV with respect to each other. As the contact is much shorter than the characteristic 
inelastic length, 

( 0, e>eV/2 
f(e,x) = i 1-x/L, eV/2> e> -eV/2 (17) 
[ 1, e < -eV/2. 

With this distribution function, the expression for the spectral density of extraneous currents (0) takes the form 



{5f*\v)5]f\T>))^ = 4aS a0 5(r - r')| (l - f ) 



(18) 
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III. CIRCULAR-SECTION CONTACT. ANALYTICAL RESULTS 



Consider the Poisson equation with the boundary conditions (|13|). As the system is axially symmetric, all the 
quantities may be considered as independent of the azimuthal angle and dependent only on the longitudinal coordinate 
x and radius r. In this case, the boundary condition (|l3|) takes the form 



r °~dr~ f 



= 0, 



-uo e d r 
-iui + 47TCT 8q 



(19) 



Supposefirst that p is real and positive. Then one may introduce a system of normalized eigenfunctions tj) n satisfying 
the equation 



with the boundary conditions (|l^). These functions are given by 

1 Jo(k n r) 



v 1,2 r y/J^(k n r ) + Jf(k n r )' 



(20) 



(21) 



where Jo and J\ are the Bessel functions of zeroth and hrst order and the eigenvalues k n are determined from the 
equation 



, J\{k n r Q ) 

knro-r-r, r = M- 

Jo{k n roj 



(22) 



Since functions ip n form an orthogonal basis, for an arbitrary charge-density fluctuation Sp, the solution of Poisson 
equation (|To|) with the boundary condition (|l9|) is given by the expression 



oo „ „ 

x,r) = -4:ir^2ip n (x,r) / dx g n (x, x') dS' ip n (r')p(x' ,/), 
n=0 { J 



(23) 



where dS' = 2itr'dr' and g n (x,x') is the Green's function of the equation 

with the boundary conditions g n (0,x') = g n {L. x') — 0, which is given by the expression 

sinh(fc„a;) sinh[fe n (L — x')] 



(24) 



g n (x < x') 
g n (x > x') = - 



k n sinh(fc I jL) 
sinh(fc„a;') sinh[fc„(L — x)] 
k n sinh(fc„L) 



(25) 



Equation ( p3|) may be analytically continued to complex values of p given by Eqn. ( |l9| ) . Substituting Eqn. (Q) for 
8p into Eqn. (p3[) for 8<p and then ( p3| ) into (|T6|), one obtains the expression for the fluctuation of the current flowing 
through the left end of the contact in the form 



5/(0) = S ^>2i^ f dx' f dS' 

n=0 I J 



d 2 g n (x 1 x') 



dxdx' 



dg n (x,x') 



x=0 



Ox 



x=0 



Or* Jr 



(26) 



where Sq = nr^ and tp n is tp n averaged over the cross section of the contact: 



"~>0 



2J 1 (k n r ) 



■K l / 2 rlk n ^J§(k n r ) + Jl(k n r 



(27) 
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To obtain the fluctuation of the current flowing through the right end of the contact, SI(L), one must substitute 
x = L for x = in Eqn. (p6|). 

At oj = 0, all transverse modes with n^O have vanishing cross-sectional averages, tp n = 0, and the corresponding 
longitudinal factors g n exponentially decay at \x — x'\ > r$. This is quite natural because the electrical field produced 
by a charge inside the contact cannot penetrte through its outer surface and is uniformly distributed over the contact 
cross section at large distances from the source. Hence the only contribution to Eqn. (|26|) will be given by the lowest 
transverse mode with ko = and iJjo(t) — n^ 1 / 2 ^ 1 . In this case, Eqn. ( p6f ) takes the form 

L 

61(0) = LJdxJ dS8f x xt . (28) 
o 

This is just the result obtained in ||. 

Consider now the case where the contact length L is much larger than its diameter 2ro and the frequencies are 
sufficiently low, i.e., to <C AiraSo/sdro- In this case, the corrections to the zero-frequency eigenfunctions ip n , as well as 
the corrections to the products k n ro with n ^ 0, are proportional to /i and therefore small; hence the contributions 
to 81 ( |2^ ) from the modes with n^O remain insignificant. However, the lowest eigenvalue is given by 

fc = r - 1 (2 M ) 1 / 2 , (29) 

and the product k$L may be sufficiently large. Therefore, the contribution from the lowest mode governed by ga(0, x) 
may change significantly. In view of this, the expression for 81 takes the form 

8 m =J dxko COSh ^- X)] [USX?, (30) 
J smh(k L) J 

o 

where ko is given by Eqn. (p9|). Physically, this implies that the contact is represented as an alternating series of 
resistors with generators of random current and grounding capacities connecting the electrodes (see Fig. ||). Note that 
81(0) is phase shifted with respect to the extraneous current inducing it. Multiplying Eqn. (pQ) by its complex con- 
jugate, substituting the spectral density of extraneous currents ( |l8| ) into the product, and performing the integration 
with respect to x, one obtains: 



Sf L (w) = 2eJ 

where 



1 sinh(27 w L) - sm(2-/ UJ L) 
jujL cosh(27„L) - cos(27 w L) 



(31) 



* = y^Rr- (32) 

The frequency dependence of the shot noise is shown in Fig. ||. At zero frequency, we rederive the well known result 
Sj L = |e/. However, at frequencies about aSoro/sdL 2 , the spectral density sharply rises and tends to the full value 
of classical shot noise, Sj L — 2el. This suggests that the corresponding correlation function is negative at t ^ t'. The 
anticorrelation is the consequence of the Coulomb repulsion of electrons: an entrance of an electron into the contact 
decreases for some time the probability for another electron to enter it, similarly to the case of a single-electron 
transistor @, |J. 

Along with the spectral density of noise at one end of the contact, one may also consider the cross-correlated 
spectral density 

Sf R = \ (81(0, u)5I(L, -u) + 81(0, -u)8I(L, u)) , 
which describes the correlation between the currents flowing through the opposite ends of the contact. Multiplying 



Eqn. (30) by its complex conjugate for 81 (L) and performing the integration with the spectral density of extraneous 



currents (Uq), one obtains: 



S LR ( ) - cosn (^-k) sin(7^L) - cos(7^L) sinh^L) 
1 [UJ) ~ lu L cosh(2 7w L) - cos(2 7u L) ' [ ' 
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The frequency dependence of Sf R is also shown in Fig. ||. At w = 0, it also equals tel. However in contrast to 
Sj L (uj), it sharply decreases with increasing frequency and tends to zero in an oscillatory way with further increase 
of frequency. 

Consider now the high-frequency limit. In this case, the boundary condition ( fl3| ) takes the form ?/vi( r o) = 0, so 
that the quantities k n ro are the zeros of zero-order Bessel function. In this case, functions ?p n are real and form an 
orthogonal system. Owing to the orthonormality conditions, the expression for the spectral density of noise may be 
written in the form 



L 

Sf L (oo) = AS 2 eVaJ2^K 2 j dx' |(l - ^ 



{ 


d 2 g n (x,x') 




2 

+ k 2 n 


dg n (x,x') 




') 




dxdx' 


x=0 


dx 


x=0 





(34) 



As Joik n r ) = 0, Eqn. ( |27| ) reduces to ip n = 2ir 1//2 r 2 k n 1 . Substituting the explicit expressions for g n ( p5| ) into 
Eqn. (34), one obtains: 



°° 1 

Sf L (oo) = &elJ2jT sa coth(fc„L) 

n=l ^ nr °> 



Similarly, one obtains for the cross-correlated spectral density: 

1 



coth(fc„L) 



coth(fc„L) 



k n L 



^ (k n ro) 2 sinh(k n L) ^"""V"-"; 
In the limiting case of ro 3> £, both expressions take the form 

In the opposite limiting case of r <C Eqn. ( |35| ) takes the form 



(35) 



(36) 



(37) 



(38) 



whereas Sf^oo) ( J3q ) tends to zero according to the exponential law. The L/tq dependences of Sf L (oo) and Sf R (oo) 
are shown in Fig. HP 



IV. PLANAR CONTACT. NUMERICAL RESULTS 



Consider now a planar contact in the shape of a layer of thickness do i n the y direction (0 < y < do) and of width 
W (W 3> max(d ,L)) in the z direction, the average current flowing in the x direction. Because of large W, the 
effects of boundaries in the z direction may be neglected and all the quatities may be considered as independent of z. 
Introduce an orthonormal system of functions 



sin(g n x), 



9n = 



(39) 



which obey the boundary conditions y n (0) = (p n (L) = 0. For an arbitrary charge-density fluctuation dp, the potential 
fluctuation 6cj> induced by it may be presented in the form 



do 



x,y) = -4n^y n {x) / dx' tp(x') / dy' Q n (y,y')5p{x' ,y') 



(40) 



G 



where Q n (y,y') satisfies the equation 



with the boundary conditions 
/ -iu)£ d 



Qn 



dQ n 
dy 



= 0, 



y=d 



A-koSq n dy 



= 0. 



(41) 



(42) 



y=o 



Explicitly, Q n for y > y' is given by the expression 



, n _ 1_ q n do cosh[g„(d - y)} - z»sinh[g„(rf - y)} 

2q„ <j„d cosh(<? n d /2) - iQ,smh(q n d /2) 

q n d cosh(q„y') - ifl sinh(g„y') 



q n d sinh(g„d /2) - cosh(q„d /2) ' 



(43) 



where — u>Sddo/4TtaSo is the dimensionless frequency. The corresponding expression for y < y' is obtained from 
Eqn. ( |4^ ) by interchanging y and t/'. Substituting Eqn. @ for (5p into Eqn. J40| ) and then substituting ( |40| ) into (|lt{), 
one obtains the expression for the fluctuation of current flowing through the left end of the contact in the form 



d d 



61(0) = WJ2^T fdx'fdyfdy' 

„=! dX *=0{ { { 



dy' 



-&3l 



(44) 



The fluctuation of current flowing through the right end of the contact may be obtained by substituting x = L for 
x = in Eqn. (|44|). Using the correlator of extraneous currents (|l8|), one obtains the expressions for the spectral 
densities Sj L and Sj R in the form 



Sf L = 8eld J2 kmkn ( M "« P ™« + M 'Jan P mn) , 
m— 1 n— 1 

OO OG 

S\ R = ieldo ]T EK- 1 )" 1 + (-l) n ]*m*» (AW^, + M^ n P mn ) 



(45) 
(46) 



m— 1 n— 1 

In these expressions, we used the notation 



do do 



M mn = -^ J dy J dyx J dy 2 Q m (yi,y)Q n (y2,y), 



(47) 



do do do 



M 'n\n ^ h j ^ j j ^ 







dQ m (yi,y) dQ* n {y 2 ,y) 
dy dy 



(48) 



dx(p m {x)ip n {x)^ (l - ~J , (49) 



dip m dip n X 

dx dx 



(50) 
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Using the notation t m — tanh(fc m <io/2) and D m = (k^d^t^ + fl 2 ) 1 and performing partial summation over the 
internal index, one may bring Eqns. ( |45|) and (|46| ) to the form 

SPM _ a rf + ^ f { (| _ ,i^r) ^ - - 4) 

oo / D D 1 

Fo' / [1 ~l~ ( — 1) ~^ ] j 7 ; o j o\9 {dnk m t m kntn + fi ){k m tm + k n t n ) >, (51) 



8 °° ' D D 1 

Fo / [1 "F (— l)'" - '"™] , „ „>„ {dnk m t m k n t n + f2 2 )(fc m t m + k n t n ) > 7 (52) 

where the primes by the sums over n show that n ^ m. 

The contour plots of 5*f L and Sf^ vs. logarithms of frequency and contact length are shown in Figs. || and ||. 
Qualitatively, their behavior is similar to that in the case of a circular-section contact. At low frequencies and small 
contact lengths, both quantities tend to |e/. At high frequencies and large contact lengths, Sj L and Sj R tend to 2el 
and zero, respectively. It is also clearly seen that at L/do > 2.86, the frequency dependences of Sf R exhibit negative 
portions. 



V. CONCLUSION 



Both circular and planar contacts exhibit qualitatively similar noise properties. At small length-to- width ratios, 
when the screening of charge fluctuations by the electrodes is more efficient than the screening by the ambient medium 
and pile-up of the charge in the contact is forbidden, the effects of long-range Coulomb interaction reduce to averaging 
the extraneous currents over the volume of the contact at arbitrary frequencies. The situation is different, however, 
for long and narrow contacts, where the charge fluctuations are mostly screened by the ambient medium and pile-up of 
charge in the contact is allowed. At sufficiently high frequencies, the correlation length of fluctuations becomes smaller 
than the length of the contact. In this case, the fluctuations of current at the ends of the contact, which are observed 
in the external circuit, are dominated by extraneous currents in the narrow adjacent layers. The corresponding 
spectral densities are equal to that of the classical shot noise, 2el, while the fluctuations at different contact ends are 
completely independent. 

This work was supported by DOE's Grant #DE-FG02-95ER14575 and by the Russian Foundation for Basic Re- 
search (project #96-02-16663-a). 
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FIGURE CAPTIONS 
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FIG. 1. Longitudinal cross section of the contact. The dotted rectangle is the metal with impurities, thin solid lines show 
the contact - electrode interfaces, thick lines show the dielectric layers of thickness So, and the hatched areas show the grounded 
ambient medium. 

FIG. 2. Physical model of noise in a long and narrow contact. Each section of the R — C line contains a generator of random 
current. 

FIG. 3. Dependences of the normalized spectral density of noise at one of the contact ends Sf L /2eI (solid line) and the 
cross-correlated spectral density Sf R /2eI (dashed line) on the dimensionless frequency uiL 2 Ed/4iv8oro for a long narrow contact. 

FIG. 4. Dependences of the normalized spectral density of noise at one of the contact ends Sf L /2eI (solid line) and 
cross-correlated spectral density Sf R /2eI (dashed line) on the length-to-radius ratio L/ro in the high-frequency limit. 

FIG. 5. Contour plots of Sf L vs. logarithms of normalized frequency Q = LO£ d do/4:ira6o and normalized length L/do for 
the planar contact. 



FIG. 6. Contour plots of Sf R vs. logarithms of normalized frequency Q — cj£ddo/4na5o and normalized length L/do for 
the planar contact. 



9 



